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We assess the current phenomenological status of transverse momentum dependent (TMD) parton 
distribution functions (PDFs) and fragmentation functions (FFs) and study the effect of consistently 
including perturbative QCD (pQCD) evolution. Our goal is to initiate the process of establishing 
reliable, QCD-evolved parametrizations for the TMD PDFs and TMD FFs that can be used both 
to test TMD-factorization and to search for evidence of the breakdown of TMD-factorization that 
is expected for certain processes. In this article, we focus on spin-independent processes because 
they provide the simplest illustration of the basic steps and can already be used in direct tests 
of TMD-factorization. Our calculations are based on the CoUins-Soper-Sterman (CSS) formalism, 
supplemented by recent theoretical developments which have clarified the precise definitions of the 
TMD PDFs and TMD FFs needed for a vahd TMD-factorization theorem. Starting with these 
definitions, we numerically generate evolved TMD PDFs and TMD FFs using as input existing 
parametrizations for the coUinear PDFs, coUinear FFs, non-perturbative factors in the CSS factor- 
ization formalism, and recent fixed-scale fits. We confirm that evolution has important consequences, 
both qualitatively and quantitatively, and argue that it should be included in future phenomeno- 
logical studies of TMD functions. Our analysis is also suggestive of extensions to processes that 
involve spin-dependent functions such as the Boer-Mulders, Sivers, or Collins functions, which we 
intend to pursue in future publications. At our website [l|] we have made available the tables and 
calculations needed to obtain the TMD parametrizations presented herein. 



I. INTRODUCTION 

The factorization theorems of pQCD have been in- 
strumental in the successful application of QCD theory 
to phenomenology. The standard coUinear factorization 
formalism makes use of "integrated" PDFs and FFs 
which depend only on a single longitudinal momentum 
fraction, while the small momentum components, includ- 
ing the transverse components, are integrated over in the 
definitions. The integrated PDFs and FFs have consis- 
tent operator definitions in QCD, with appealing inter- 
pretations in terms of parton model concepts. However, 
the standard coUinear factorization formalism relies on 
approximations that are only valid for sufficiently inclu- 
sive observables. In order to address many of the issues 
now at the forefront of research in QCD and its role in 
hadron structure, pQCD factorization must be extended 
to situations where the usual approximations are not ap- 
propriate. 

A transverse momentum dependent TMD- 
factorization formalism goes beyond the standard 
factorization framework by allowing the PDFs and 
FFs to depend on intrinsic transverse momentum in 
addition to the usual momentum fraction variables. 
As such, different sets of approximations are needed 
in the factorization proofs. The PDFs and FFs in a 
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TMD-factorization formalism are referred to as TMD 
PDFs and TMD FFs (they are also called "unintegrated" 
or " /cj^-dependent" ) to distinguish them from the more 
familiar integrated correlation functions of coUinear 
factorization. Henceforth, we will refer to TMD PDFs 
and TMD FFs coUectively as "TMDs" . 

The role of the intrinsic transverse momentum car- 
ried by partons in high energy collisions is becoming 
increasingly central in discussions of how to probe the 
details of hadronic structure in high energy collisions. 
The usefulness of the TMD concept is evident from the 
large variety of situations where it makes an appear- 
ance. Generally, TMD-factorization is needed to describe 
processes that are sensitive to intrinsic parton trans- 
verse momentum. The Drell-Yan (DY) process, single- 
inclusive deep inelastic scattering (SIDIS), and back-to- 
back hadron production in electron-positron annihilation 
at small transverse momentum are all classic examples of 
where TMD-factorization formulae are frequently used. 
More recently, TMD PDFs and FFs have been under 
intense study as objects that carry information about 
the spin structure of hadrons; the Boer-Mulders, Sivers, 
and pretzelosity functions are all specific examples of 
TMD PDFs, while the Collins function is an example of 
a TMD FF. For a recent review of TMDs in spin physics, 
see Ref. [3]. In very high energy (small-x) resummation 
physics, where there is a lack of /cT-ordering, the TMD 
gluon distribution is especially important, and similar 
issues must be dealt with. Finally, TMD functions are 
useful tools in the construction of Monte Carlo event gen- 
erators, where the details of final state kinematics are of 
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interest. 

While TMDs can potentially provide a much deeper 
understanding of QCD and hadron structure, the the- 
oretical framework of TMD-factorization is much more 
complicated than the more standard collinear factor- 
ization. In derivations of collinear factorization, there 
are important cancellations that occur after integra- 
tions of parton momentum are carried out. With TMD- 
factorization, the integrals over parton transverse mo- 
mentum arc left undone in the definitions of the TMDs, 
and contributions that would ordinarily cancel in a 
collinear factorization treatment must be accounted for. 
Colhns, Soper and Sterman (CSS) constructed a TMD- 
factorization formalism that deals with the main 
complications of transverse momentum dependence, and 
provides a systematic treatment of pQCD over the full 
range of transverse momentum. The CSS formalism has 
proven highly successful in specific phenomenological ap- 
plications such as in the calculation of the transverse 
momentum distributions in DY processes (see, for exam- 
ple, (3,0), and is also well-suited for the production of 
back-to-back particles in e^e~ annihilation. The same 
methods are needed for the discovery of new particles 
like the standard model Higgs boson [9l-[TT|. Further- 
more, extensions of the CSS TMD-factorization formal- 
ism have been derived for SIDIS [T2l - [l^ . and including 
spin in Ref. [isj . 

However, the most common methods for applying the 
CSS formalism are not ideally suited for studies that 
are specifically oriented toward understanding the TMD 
PDFs and TMD FFs themselves. Furthermore, the re- 
lationship between the full pQCD treatment of factor- 
ization and parton-model intuition has remained much 
less clear in TMD-factorization than in collinear factor- 
ization. This has led to considerable confusion about 
how the study of TMDs should be approached in pQCD. 
That confusion is especially apparent from a comparison 
between current applications of TMD-factorization and 
collinear factorization: While there have been extensive 
programs dedicated to parametrizing and evolving the in- 
tegrated PDFs and FFs (making them indispensable and 
portable tools for phenomenology), a generally agreed 
upon framework for dealing with TMDs in an analogous 
way has not yet been established. More disturbingly, 
there has been a persistent lack of clarity or agreement 
regarding the definitions of the TMDs. A suitable set of 
definitions must be dictated by the requirements of fac- 
torization and universality, but the most common and 
naive definitions lead to inconsistencies, including un- 
physical divergences. 

Over roughly the past decade, there has been steady 
progre ss toward a better understanding of what is needed 
p^| - |25| . The issue of finding the right definitions has now 
been es pecially brought into focus by the recent work of 
Collins Mug. The definitions for the TMDs in Ref. ^ 
are uniquely determined by the requirements of factoriza- 
tion, maximal universality, and internal consistency. (By 
"maximal universality" we mean that the same correla- 



tion functions appear in a large number of processes.) 
The confusion over definitions therefore appears to be 
solved. With the new definitions, the implementation of 
evolution via the CSS formalism is not modified signifi- 
cantly from earlier treatments which means that existing 
parametrizations of the non-perturbative parts can still 
be used. However, they can now be understood as contri- 
butions to separate, QCD-evolved TMDs. Moreover, as 
we will discuss in the next two sections, the new defini- 
tions have a much more direct relationship with more 
intuitive, parton model based ways of viewing TMD- 
factorization. 

There has been much recent work devoted to 
parametrizing TMDs by assuming a parton model pic- 
ture of TMD-factorization and directly fitting cross sec- 
tion calculations to experimental data [28.~33j . This ap- 
proach to TMD phenomenology is often called the gen- 
eralized parton model (GPM) [34| . 

In addition, by working with non-perturbative models, 
it is possible to study the general properties of the TMDs 
and their relationships with each other. (See Ref. fsl] 
and references therein for an overview of this subject.) A 
famous example is the illustration via model calculation 
that the Sivers function is non-vanishing in SIDIS (sgI . 

However, most efforts to parametrize or model TMDs, 
particularly the spin-dependent TMDs, have ignored the 
role of evolution. For the unpolarized TMD PDFs 
there have been detailed implementations of evolution 
(e.g. [ll,!!^) However, even in the unpolarized case, the 
identification of the separate evolved TMDs and their 
relationship with the fundamental definitions in a com- 
plete treatment of factorization has remained unclear. 
The main purpose of the present article is to initiate 
the process of consistently including QCD evolution in 
parametrizations and models of TMDs by following the 
definitions in Ref. [l^. Addressing evolution is now an 
especially pressing task, given the very wide range of en- 
ergy scales set to be probed in experiments in the near 
future, from Jefferson Lab (JLab) to the LHC. It has 
already been shown in Refs. [i^, |4l| that evolution (in 
the form of Sudakov suppression) should be expected to 
be large. In our analysis, we will illustrate how this fol- 
lows from the evolution of the individual TMDs. We also 
aim to facilitate the future implementation of evolution 
in studies of TMDs by clarifying the relationship between 
the parton model description of TMD-factorization and 
the CSS formalism. To this end, we have made computa- 
tions available at [l[ that illustrate how to obtain evolved 
TMDs given a choice of non-perturbative starting input. 

An alternative approach to probing spin effect is to 
consider higher twist collinear functions such as the Qiu- 
Sterman function (42j . By taking transverse momentum 
moments of cross sections it is possible to relate TMDs 
to these higher twist functions [l^. See Ref. Q for 
recent work on the evolution of weighted spin-dependent 
correlation functions. 

Establishing reliable, evolved TMDs is also important 
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for testing factorization and searching for instances of 
factorization breaking. Discussions of non-universality 
or factorization breaking often bring in the concept of 
Wilson lines (or gauge-links), and the process depen- 
dent properties that can be associated with them. Al- 
ready when comparing the Sivers function in SIDIS and 
the DY process, one must account for a well-known sign 
flip that is due to the reversed direction of the Wilson 
lines in these two processes [l^]. More recently, it has 
been shown that obtaining consistent Wilson lines is even 
more problematic in the hadro-production of hadrons or 
jets {Hi + H2 ^ H3 + H4 + X). In such cases that 
require TMD-factorization, it was found that at a mini- 
mum the Wilson lines for se par ate TMDs are highly com- 
plex and process dependent |45l - l47l | . It was later shown in 
Ref. |48| that TMD-factorization generally breaks down 
completely in the hadro-production of hadrons. That is, 
separate TMDs cannot even be defined for each external 
hadron regardless of what Wilson lines are used in the 
definitions. 

The complication in the case of the complete break- 
down of factorization is caused by a failure of the 
usual gauge invariance/Ward identity arguments that are 
needed in a factorization proof. A confirmation of this 
effect would point to interesting new features of strong 
interaction physics, given that the breakdown of factor- 
ization is in just the range of kinematics where factoriza- 
tion would naively be exp ected to hold. Calculations in a 
GPM framework [49| - |51| will be needed for making pre- 
dictions that can be compared with experiment to test 
factorization and/or search for factorization breaking. 
Furthermore, computations using the methods in, for ex- 
ample, Refs. [52- 54j can potentially help to quantify and 
better understand the factorization breaking mechanism. 
It may soon be possible to find experimental evidence 
for factorization breaking, particularly in the analysis of 
RHIC data. (See, for example, the recent analysis of 
Ref. [55l - l58| .) However, definitive tests of factorization 
or factorization breaking can only become possible with 
a more precise determination of the TMDs. Another mo- 
tivation for this article is therefore to begin the deter- 
mination of TMD parametrizations that should be used 
in tests of factorization in the case of spin-independent 
hadro-production of back-to-back hadrons. Finally, in 
this paper we include some details of the calculation of 
the evolution of the TMD PDFs that did not appear in 
Ref. dl 

The paper is organized as follows: in Sect. |lTl a brief 
background of TMD-factorization is provided, and the 
main complications that arise in the context of pQCD 
are discussed. In Sect. lIIll we set up the notation, and in 
Sect. IIVI we explain the definitions of TMD correlation 
functions. In Sect. |Vl we discuss the evolution of the 
TMDs in terms of the CSS formalism. We apply evolu- 
tion to existing unpolarized quark TMD fits in Sect. IVI[ 
and we present some numerical results. We conclude with 
an overview and a discussion of future work in Sect. I VIII 
In the appendices, we provide some details of the pertur- 



bative parts of our calculations. 



II. THE DIFFERENT PICTURES OF 
TMD-FACTORIZATION 

In this section we expand on the general remarks in 
the introduction by providing a very schematic overview 
of the different ways TMD-factorization is approached 
in phenomenological applications. We discuss the rela- 
tionship between parton model intuition and full pQCD 
while emphasizing the complications that can arise. We 
explain the potential for confusion when evolution and 
soft factors are included, and how this is solved by the 
use of fully consistent definitions for the TMDs. 



A. The Generalized Parton Model 

We start by considering the simplest parton model pic- 
ture of high energy collisions. There, the concept of a 
TMD-factorization formula becomes very intuitive and 
easy to state. The cross section is simply a partonic 
sub-process, folded with TMD PDFs and TMD FFs. In 
SIDIS, for example, the hadronic tensor is written as 

/ 

J cfkiT <fk2T Ff/p{x, kir) Di-,/f{z, zk2T)^ 

X(5(2)(kiT + qT-k2T)- (1) 

Here \Hf{Q)'^\'^'^ describes the hard partonic sub-process, 
7*5 — q, for scattering off a quark of flavor / as a func- 
tion of the hard scale Q. (It also includes any overall fac- 
tors needed to make the left side a proper hadronic ten- 
sor.) The size of q-r is a measure of the non-collinearity 
in the process. Within the parton model, the TMDs 
Ff/p{x, kir) and Df^^f^z, zk2T) have simple probabilistic 
interpretations; Ff /p{x,kiT), for example, is the proba- 
bility density for finding a quark of flavor / with momen- 
tum fraction x and transverse momentum ki^ inside the 
proton. Equation (P) is closely analogous to the standard 
collinear factorization theorem of inclusive processes [59{ . 
The only difference is that the TMD PDFs and FFs are 
allowed to carry intrinsic transverse momentum. 

The intuitive approach to TMD-factorization embod- 
ied by Eq. ([1]) forms the basis of much current TMD phe- 
nomenology. However, derivations of TMD-factorization 
order-by-order in real pQCD involve complications that 
are not immediately apparent in parton model reasoning. 
Indeed, some aspects of partonic intuition can quickly 
lead to incorrect results if taken too literally. Further- 
more, without a complete derivation of factorization, the 
issue of universality (or non-universality) of the TMDs in 
Eq. ID) is much less clear. 
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B. Divergences and Soft Factors 

One issue is the appearance of extra divergences. In 
addition to the standard UV divergences associated with 
the renormahzation of the theory, the TMD correlation 
functions also contain so-called "light-cone" divergences. 
They arise if the Wilson lines (gauge links) in the defini- 
tions of the TMDs point in exactly light-like directions. 
We should note that the same light-cone singularities are 
not present in the ordinary collinear correlation functions 
because they cancel in a sum over final state interactions, 
which is possible due to integrations over the loop mo- 
menta, including the kx integrals. Light-cone divergences 
correspond to gluons moving with infinite rapidity in the 
direction opposite the containing hadron, and are not reg- 
ulated by the use of an infrared cutoff, so they amount 
to a real inconsistency in the definitions of TMDs. The 
most naive definitions, therefore, are invalid, and modi- 
fications are needed for a reliable factorization theorem. 
(See Refs. [U [2^ and references therein for a review of 
this and other subtleties involved in defining consistent 
PDFs.) The light-cone divergences need to be regulated, 
typically by tilting the Wilson lines slightly away from 
the exactly light-like directions. 

Furthermore, in a TMD-factorization formula the role 
of soft gluons becomes important. (In this paper, "soft" 
refers to nearly on-shell gluons with rapidity intermediate 
between the rapidities of the colliding hadrons and pro- 
duced jets.) They imply that another correlation func- 
tion, a separate soft factor, should be inserted into fac- 
torization formulae like Eq. ([T|), in addition to the usual 
TMD PDFs and FFs. Already, the appearance of a soft 
factor seems to contradict the basic parton-model intu- 
ition. 

The complications listed above, as well as a consis- 
tent matching to collinear factorization at large trans- 
verse momentum, are accounted for in the CSS formal- 
ism 0-01 • With the Wilson lines tilted to remove light- 
cone divergences, the factorization formula acquires new 
arbitrary parameters. Predictive power is then recov- 
ered by a kind of generalization of renormahzation group 
techniques. The resulting evolution equations may be 
thought of as describing the variation of the TMDs and 
the soft factors with changes in the degree of tilt of the 
Wilson lines. Physically, this corresponds roughly to a 
variation with respect to a cutoff on the phase space avail- 
able for recoil against soft gluons. 



C. Confusion Over TMD Definitions 

While the CSS formalism has been very useful for 
past phenomenological studies, the usual implementa- 
tions bear little surface resemblance to the generalized 
partonic picture we started out with in Eq. ([1]). For ex- 
ample, in Ref. (soj (and similar applications of the CSS 
formalism to the DY process) the effects of evolution are 
gathered into separate factors, and it is not clear how 



they relate to separate TMD PDFs. In other treatments 
(e.g. [13, [lO]), factorization formulae for SIDIS are pro- 
vided which contain explicit evolved TMDs, but they also 
involve separate explicit soft factors, and the hard part 
has explicit dependence on light-cone divergence cutoff 
parameters. Moreover, given the general observations 
that are reviewed in Refs. [2l|, [2^ |6l| , it is questionable 
whether the most commonly quoted definitions of the 
TMDs are even fully consistent. 

The original work of Collins and Soper Q used a non- 
light-like axial gauge to regulate the rapidity divergences. 
Later, Collins and Hautmann [l^, [13, i21,] proposed def- 
initions in which the main legs of the Wilson lines are 
light-like, but in which there is a division by a special 
type of soft factor which cancels the rapidity divergences. 
However, these definitions continue to suffer from prob- 
lems, including the appearance of badly divergent Wil- 
son line self-energy contributions as discussed recently in 
Ref. [25^ and also utilized in the treatment of TMD PDFs 
in Ref. ^24]. While many of these issues have typically 
been discussed in the context of TMD PDFs, the same 
problems arise in the treatment of FFs. 

Finally, it has remained unclear how the TMDs that 
have been used in past applications of evolution and the 
CSS formalism are related to the TMDs of other ap- 
proaches, such as those based more on generalized par- 
ton model pictures. In parametrizations of TMDs, the 
role of the soft factor is often not explicitly included 
and evolution is ignored. Many other theoretical TMD 
studies continue to quote definitions with exactly light- 
like Wilson lines. Knowledge of the operator definitions 
for the TMDs is also needed for lattice TMD calcula- 
tions [H, [g^, and in model calculations. Clearly, a more 
unified treatment of TMD-factorization is necessary in 
order to bring together these different approaches to the 
study of TMDs. 



D. Consistent Definitions, TMD-factorization and 
Evolution 



What is needed, in addition to fully consistent TMD 
definitions, is a formulation that retains as much as pos- 
sible the basic factorized structure of Eq. ([T]) , but which 
appropriately includes evolution and the effects of soft 
factors. Ideally, the situation should be closely analogous 
to what already exists for collinear factorization. Namely, 
it should be possible to clearly identify consistent and 
universal TMDs that can be tabulated or parametrized 
and then reused in formulae like Eq. ^ to make predic- 
tions for a wide variety of processes. Fortunately, this 
has been achieved in the recent work of Ref. [26| . There, 
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the factorization formula for SIDIS takes the form: 

/ 

J (fkiT (fk2T Ff/p{x, kiT; Cf) Dh/f{z, zk2T; ^J■^, Cd)x 

X d^^\kiT + qr - k2T) 

+ Y{Q,qT)+0{{A/Qr). (2) 

The first term on the right-hand side of this equation 
(responsible for the low-q^ behavior) has exactly the 
structure of the partonic TMD-factorization formula in 
Eq. ([T|), apart from the scale dependence denoted by 
/i, Ci? and Cd- The arguments Cf and Cd will be dis- 
cussed more in the explanation of the TMD definitions 
in Sect. IIVI They are left over from the need to regu- 
late light-cone divergences, and should obey ^/CfCd ~ 
0{Q^). In terms of more familiar variables, they are de- 
fined as: 



and 



(3) 



(4) 



Here, x and z are the usual Bjorken scaling and frag- 
mentation variables, Mp is the proton mass and Mh is 
the mass of the produced hadron. The rapidities of the 
proton and produced hadron are yp and yh respectively. 
The rapidity j/s is an arbitrary low-rapidity cutoff param- 
eter that separates partons with large forward rapidity 
(in the proton direction) from backward rapidity (in the 
produced hadron direction) . Variations of these functions 
with j/s will be determined by the evolution equations. 

The scale /i is the standard renormalization 
group (RG) scale. The TMD correlation functions, 
f//p(x, kiT;Ai;C-F) and ^/^//(z, zk2T; m; C-d)i have 
definite and consistent operator definitions. They 
include the effects from soft gluons in such a way that 
no soft factor appears explicitly in Eq. Evolu- 
tion can be implemented on Fj/p(a;, ki^; Cf) and 
Dh/ zk2T', H'iCd) independently, and the basic steps 
closely follow the usual CSS approach. We will discuss 
the definitions more in the next section, but for now we 
mention that they solve most of the theoretical problems 
summarized in Refs. [2l|, and Sect. IlICi including 
the appearance of light-cone divergences and Wilson line 
self- interactions. 

The term, F((5,gT), accounts for the large-q^ de- 
pendence of the cross section, where the approxima- 
tions needed for TMD-factorization break down. There, 
collinear factorization becomes the appropriate frame- 
work. The error term is suppressed by [K/QY where 
a > 0. The first term on the right side of Eq. ([2]) is valid 
up to corrections of order {qt/Q)"', but the Y{Q,qT) is 
needed for a valid treatment of factorization over the full 
range of qy. 



The derivation of Eq. ([2]) within pQCD factoriza- 
tion, with consistent definitions for the TMDs, is an im- 
portant breakthrough because it connects TMD studies 
from a GPM framework with formal QCD and clarifies 
the meaning of TMD evolution. We will use Eq. ([2]), 
alon g wi th the associated definitions for the TMDs from 
Ref. [26], to obtain momentum space fits for use in phe- 
nomenology. The non-perturbative input can be ob- 
tained from already existing models or fits made at fixed 
scales. For the TMD PDFs, much information about the 
non-perturbative input is already available from fits that 
use the standard b^-space formulation of the CSS for- 
malism in the DY process. 



III. SETUP AND NOTATION 

We start by setting up the basic notation. In our 
convention for light-cone variables, a four-vector = 
{V'^ ,Yt) has components. 



Vt - {v'',vy). 



(5) 



The z-component picks out the forward direction. Note 
that ^2 = 2V+V- - Vf,. 

For the processes we are interested in, there are always 
two relevant light-like directions which we label ua and 
ub, and define to be: 



MA = (l,0,Ot) UB = (0,l,Ot). 



(6) 



In the SIDIS example, ua and ub characterize the direc- 
tions of the incoming proton and the produced jet. A 
Wilson line from a coordinate a; to oo along the direction 
of a four- vector n is defined as usual: 



W{oo, x;n) = P exp 



-Wo 



dsn- Aq{x + 371)1" 



(7) 

In these definitions, the bare fields and couplings are 
used, P is a path-ordering operator, and t°- is the gener- 
ator for the gauge group in the fundamental representa- 
tion, with color index a. 

As discussed in the previous section, light-cone diver- 
gences must be regulated by tilting the direction of the 
Wilson line away from the exactly light-like direction. 
Therefore, we need to define another set of vectors 
and riB analogous to Eq. ^ but slightly tilted, so that 
they have rapidities yA and ys- 



= (-e'^^l,0^ 



(8) 



Note that the tilted Wilson line directions are space-like, 
n\ ^ nl < 0. The use of space-like directions for the 
Wilson lines ensures maximum universality for the def- 
initions of the TMDs, as explained in Ref. [63|. In all 
of our calculations, fi is the standard MS mass scale in 
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dimensional regularization and the dimensional regular- 
ization parameter e is defined in the standard way as 
2e — A — d where d is the dimension of space-time. 

Though our results apply generally to the standard 
factorizable processes, we will continue to use SIDIS as 
a reference process for explaining the definitions. Let us 
also rewrite the TMD-factorization formula for SIDIS in 
Eq. dD) as 

W'"" = 

Y,\■HfiQ;^lf^ [ d'kiTd'k2T^(')(kiT + qT-k2T) 
X Ff/p{x, kiT; n; Cf) Dh/f{z, zk.2T; Cr>) 

X Ff/p{x,hT;fi;CF)Dh/f{z,hT;n;CD)- (9) 

Throughout this paper, it will be implicit that all momen- 
tum components are in the hadron frame. (The hadron 
frame is where both hadrons have zero transverse mo- 
mentum and is a natural frame for setting up the steps 
for factorization.) 

Hereafter, the Y{Q, q^) term that appeared in Eq. ([2]) 
will also be dropped because our primary interest is in 
the qt << Q regime where TMD-factorization is appro- 
priate. Also, we will drop any explicit +0{{A/Q)'^) sym- 
bols. We have written the TMD-factorization formula in 
coordinate space in the second equation of © because 
it is simpler to explain the coordinate space definitions 
of the TMDs and their evolution. Later we will Fourier 



transform the TMDs back to momentum space when we 
analyze them numerically. 



IV. DEFINITIONS OF THE TMDS 

As explained in Sect. |TT1 our calculations are based on 
the formulation of TMD-factorization explained in detail 
in Ref. [26]. A repeat of the derivation is beyond the 
scope of this paper. However, in order to put our later 
calculations into their proper context, we will give an 
overview of the basic features of the formalism in this 
and the next section. We refer the reader directly to 
Ref. [1^ for pertinent details. 



A. Soft Factor Deflnition 

We have already stressed in Sect. Ill CI that the defini- 
tions of the TMDs in Eq. ^ are not the often quoted ma- 
trix elements of the form ^ {P\tp Wilson Line tp\P) with 
simple light-like Wilson lines connecting the field oper- 
ators. Using such definitions in a factorization formula 
leads to inconsistencies, including unregulated light-cone 
divergences. Also, soft gluons with rapidity intermedi- 
ate between the two nearly light-like directions need to 
be accounted for in the form of soft factors. Therefore, 
before we can discuss the definitions of the TMDs that 
will ultimately be used in Eq. (j9|), we must provide the 
precise definition of the soft factor. In coordinate space 
it is an expectation value of a Wilson loop: 



^(0) (br; 2/A, Vb) = — (0| W^(bT/2, oo; nB)la W{hT/2, oo; nA)adW{-hT /2, oo; nB)bcW{-hT /2, oo; nA)\^\Q)^o s.i.- 

, m 



We have used the vectors in Eq. ([8]) to define the direc- 
tions of the Wilson lines so that, as long as ua and ys 
are finite, the Wilson lines in Eq. (|10l) are non-light-like. 
The subscripts a, 6, c and d are color triplet indices, and 
repeated indices are summed over. The "(0)" subscript 
indicates that bare fields are used. The soft factor con- 
tains Wilson line self-interaction (S.I.) divergences that 
are very badly divergent and are unrelated to the orig- 
inal unfactorized graphs. They must therefore be ex- 
cluded, and we indicate this with a subscript "No S.I.". 
We emphasize, however, that this is only a temporary 
requirement because all Wilson line self-energy contribu- 
tions will cancel in the final definitions. Another poten- 
tial complication, pointed out in Refs. [19,, 20], is that 
exact gauge invariance requires the Wilson lines to be 
closed by the insertion of links at light-cone infinity in 
the transverse direction. However, the transverse seg- 



ments will not contribute in the final definitions of the 
TMDs (at least in non-singular gauges), so we do not 
show them explicitly in Eq. ([rU]) . Again, the final ar- 
rangement of soft factors will ensure a cancellation. 

Rather than appearing as a separate factor in the 
TMD-factorization formula, soft factors like Eq. ((TU)) will 
be part of the final definitions of the TMDs. Their role 
in the definitions will be essential for the internal consis- 
tency of the TMDs and their validity in a factorization 
formula like Eq. ©. 



B. TMD PDF and FF Definitions 

Now we turn to the definitions of the TMDs them- 
selves, starting with the unpolarized TMD PDF. The 
most natural first attempt at an operator definition is 
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obtained simply by direct extension of the collinear inte- tilted to avoid light-cone singularities. The operator def- 
grated parton distribution, though with the Wilson line inition is 



F;;|."^a;,bT;A*;yp-yB) 



^Ti-cTiD I ^e-"^ {P\i>f{w/2)W{w/2,^,nB)^^W{-w/2,cx,nB)i^f{-w/2)\P),,NoS.i.- (H) 

' ZTT Z 



This definition does not account for the overlap of the 
soft and collinear regions, so we refer to it as the "un- 
subtracted" TMD PDF. Here w = {0,w~,hT) and yp is 
the physical rapidity of the hadron. As usual, the struck 
quark has a longitudinal plus-component of momentum 
= xP^ . The non-light-like direction of the Wilson 
line is given by the np vector defined in Eq. ([S]) so that 
the light-cone divergences are regulated by the finite ra- 
pidity yp- As usual, the definition includes a trace over 
color. The c subscript is to indicate that only connected 



diagrams are included. Equation PT|) reduces exactly 
to the most naive definition of the TMD PDF when the 
light-like limit of ?/b — >■ — cx) is taken. Indeed, in the final 
definition we will take this limit, but then the role of the 
soft factors becomes important. 

While the definition in Eq. (jlip is intuitively appealing, 
modifications are needed in order to have a consistent 
definition that can be used in a factorization formula like 
Eq. © . The com plet e definition for a quark / in proton 
P, given in Refs. [H [13, is 



F//p(x,bT;M;CF) = F;;|."'^(a;,bT;//;yp - (-c^)). 



5(o)(bT;+oo,ys) 



\ ^(o)(bi 



- 00)5(0) (br; y^,- 



-Zp Zi 



(12) 



Here, the "oo" arguments for the rapidity variables in the unsubtracted PDF and the soft factors are meant in the 
sense of a limit. All field operators are unrenormalized, and Zp and Z2 are the PDF and field strength renormalization 
factors respectively. The soft factors on the right-hand side of Eq. (jl2p contain rapidity arguments ys ■ It is an arbitrary 
parameter which can be thought of as separating the extreme plus and minus directions. It will be convenient to 
express the dependence on j/s via C_f , defined in Eq. ^ . On the left-hand side of Eq. (fT2|) , the dependence on y^ is 
expressed via the dependence on C^p. 

Although we will not repeat the derivation that leads to Eq. ([12]), we remark that the definition is unique given 
the requirements that: a.) Factorization holds with maximal universality for the TMDs. b.) No explicit soft factor 
appears in the final factorization formula, Eq. ©. c.) Self-interactions of the Wilson lines, and attachments to gauge 
links at infinity cancel in the final definition, d.) The Collins-Soper (CS) equations are homogeneous. 

There is an analogous definition for the TMD FF. The unsubtracted version, analogous to Eq. ([TT|) . is 



^^"//•^(z, br; yA -yh)^Y. I^^TrcTrz,- / 



dw 
~2^ 



(0|7+ W{wl2, c^, nA)i^f{wl2)\h, X) 
X {h,X\^fi~w/2)Wi-w/2,oo,nAy\0)c,NoS.i.. (13) 



Now yh is the physical rapidity of the produced hadron or jet. The complete definition of the TMD FF with the soft 
factors included is: 



DH/f{z,hT;n;CD) 



yh) 



S{o){hT;ys,-oo) 



\ '5'(o)(bT;+oo,-oo)S'(o)(bT;+oo,2/5) 



Zd Z2 



(14) 



Again there is dependence on the soft rapidity ys- For rapidity scale ys via Eq. (|4|. 
the FF, the energy cutoff scale is related to the soft 
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Equations and ([14]) are the correct TMDs for the 
TMD-factorization formula in Eq. © as well as for e+e~ 
annihilation with back-to-back jets. Up to a flip in the 
direction of the Wilson line from future to past pointing, 
which is important for accounting for a sign flip in certain 
types of TMDs, the TMD PDF in ^ is also relevant 
for the Drell-Yan process. In this section we have clar- 
ified the meaning of the energy parameters and Cd, 
which were already discussed in Sects. |lTl IIIIl Together, 
Eqs. (131 SI) give VCfC-D ~ Q^- They are in principle ar- 
bitrary, and the full factorization formula in Eq. ^ is 
exactly independent of the choice of ys (and therefore 
and VCd). However, different choices are needed 
for each factor in Eq. ^ in order to optimize the conver- 
gence properties of the perturbation series. To obtain the 
TMDs appropriate for different scales, we must appeal to 
evolution equations which are the subject of Sect. [V] 



the soft factor into a product of two factors: 



C. The Role of Soft Factors 

Admittedly, the final definitions in Eqs. (fT^ and (fH)) 
appear rather complex. While detailed derivations are 
beyond the scope of this paper, it is nevertheless worth- 
while to make a few intuitive remarks about how these 
definitions arise in a treatment of factorization. A much 
more detailed treatment is found in Ref. 26] . 

For now we simplify the notation for the TMDs by 
dropping all arguments and symbols not directly related 
to Wilson line rapidities. The cross section can then be 
written (schematically) as. 



da = \n\ 



(-CX))) X D™«"'^(-hoo - yh) 



5(+oo, — oo) 



The F""™'^(yp - (-oo)) and i:)""™''(+cx) - yn) are the 
same "unsubtracted" TMDs from Eqs. ^ and (fT3| . 
They each describe the distribution of gluons in their 
relevant collinear direction, but they also both account 
for soft gluons with nearly zero rapidity. Therefore, the 
soft factor S'(-|-oo, — oo) in the denominator is needed 
to remove double counting. Since the Wilson lines in 
the unsubtracted TMDs are light-like in the plus and 
minus directions respectively, they also include rapidity 
divergences. Thus, the role of the S{+oo, —oo) is also 
to cancel these rapidity divergences. Although Eq. (fT5)) 
properly accounts for all soft and collinear regions and 
deals with the divergences, it is not factorized. Because 
of the rapidity divergences, it is immediately clear that 
punsuhj^y^ _ (_oo)) and i)""""'^(-Hoo - y,,) are not sepa- 
rately well-defined, and in the full formula they are en- 
tangled via the soft denominator. To get a factorized 
structure for Eq. (jisp , with each factor individually well- 
defined, a natural first step to try is simply to separate 



da= \H\'- 



S{+oo, — oo) 



(-00)) D 

X 



unsub 



(-I-00 - yii) 



5(-|-oo, — oo) 



(16) 

One is then tempted to identify the factors on either 
side of the "x" with the TMD PDF and the TMD FF. 
However, these definitions still contain uncanceled diver- 
gences. In each factor, the rapidity divergence in the 
numerator is not completely canceled by the square root 
rapidity divergence in the denominator, and new rapidity 
divergences are introduced by the Wilson line pointing 
in the opposite direction. So the next modification of 
Eq. (fT5|) is to write 



d<7 = \H{- 



^unsub(-y^ 



-oo)) 



S'(-|-oo, — oo) 



•y/s'( + 0O,-0o) |)u„sub(-^^ _ 



S'(+oo, — oo) 



S'(+oo, — oo) 



F 



unsub 



{yp - (-oo)) 



S'(+oO; — oo) 



\J S{ + OG,ys)S{ys,-QG) £)unsub^_^^ _ j^^J 



S{+oo,ys) S{ys,-'X}) ^/S'(+oo, -oo) 
= |Hpx 



S{+oo,ys) 



S'(+oo, — oo)S'(j/s, — oo) 



X <^i)""^"^+oo-y,0. 



5(2/5,-00) 



5(+oo, — oo)5(+oo, ys) 



(17) 



(15) 

After the first equality, we have simply multiplied and 



divided by y S'(+oo, — 00). After the second equal- 
ity, we have used the group relation S{yA,yc) 
S{yA,yB)S{yB,yc), which follows from the evolution 
equations for the soft factor — see, e.g., Ref. [1^ chapter 
10. (In fact, this expression should also include an over- 
all factor that depends on rapidity yB- But this cancels 
in Eq. (jl7p between the numerator and denominator and 
does not affect our argument.) This allows for a sepa- 
ration of the soft factors into pieces that have rapidity 
divergences only in the plus or only in the minus direc- 
tions, with any other rapidity divergences cut off by the 
arbitrary scale ys- By rearranging the soft factors, every- 
thing can then be grouped into the factors on the last two 
lines. In each factor, all spurious divergences cancel, and 
we arrive at the separately well-defined TMDs in braces. 
These correspond to the definitions in Eqs. (fT^ and (|14l) . 

To summarize, the light-like Wilson lines are needed 
in each separate "unsubtracted" TMD of Eq. ([T5]) . but 
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the contribution from gluon attachments to a Wilson line 
where the gluon has nearly the same rapidity as the Wil- 
son line does not correspond to any real physics. To 
cancel these spurious contributions to the cross section, 
there must be an equal number of both plus-pointing and 
minus-pointing light-like Wilson lines in the numerator 
and denominator, as is the case in Eq. (jlSp . Applying 
this same requirement to the separate TMDs (a TMD 
PDF and a TMD FF, in our case) leads uniquely to the 
definition in the last two lines of Eq. (IT71) and Eqs. (fT^ 
and (HH). Compare this with the situation in Ref. (65| . 
There, as in our Eq. ([T5|) . the needed cancellations occur 
in the full cross section expression, but not in the indi- 
vidual TMD factors. To get separately consistent TMDs, 
the steps summarized in Eq. (|17|) are needed. 



V. EVOLVED TMDS 

The evolution of the TMDs follows from their defini- 
tions, Eqs. We start with the evolution of the 
TMD PDF. The energy evolution is given by the CS- 
equation for Eq. 



ainF(a;,bT;/x,CF) 



K(hT■,^i) 



where the function K{hT, /i) is defined as, 

K{hT; fi)^-—\n\ 

2dys \S{hT; +00, ys) 



(18) 



(19) 



Equation p8|) follows directly from differentiating 
Eq. with respect to v^Cf and using the definition 
of K{hT', /i). Note that it is S{hT) rather than 5(o)(bT) 
that appears in Eq. (1191) . Thus it is important to account 
for the UV renormalization factors ZpZ2 in Eq. ([T^ . 

The RG equations for both F(a;, by; /i; Cf) and 
K{yiT\^ are also needed. They are. 



and 



din /i 



'iKigifJ-)) 



(20) 



-7f(5(m);Cf^'). (21) 



The functions 7_fs:(g(/i)) and 7_F(5(/i); Cf/m^) are the 
anomalous dimensions of K(hT;^J.) and F{x,hT; H,Cf) 
respectively. Using Eqs. (|18ll2ip . the energy evolution of 
7f can be derived: 



7F(5(M);CF/Ai^) =7f(5(m);1) 



i7K(.9(M))ln^. (22) 



At small-^T, Eq. ([12]) can itself be calculated within 
a coUinear factorization formalism [4]. Namely, it sep- 
arates into a perturbatively calculable hard scattering 



coeflicient and an integrated PDF, convoluted over mo- 
mentum fraction: 

Ff/p{x,hT;fJ.,CF) = 

= / -^^f/3i^/^^^T■XF,^I.,9{^l))fJ/p{x]^i) 

+ 0((Aqcd6t)'^). (23) 

The functions fj/p{x;^) are the ordinary integrated 

PDFs and the Cf/j{xlx, bx; Cf, M75(a*)) are the hard co- 
efficient functions, which are provided to first order in 
Appendix [A) The last term denotes the error, which 
grows large when br ^ ^qcd- 

At large 6t, the perturbative treatment of the bx- 
dependence is no longer reliable. In momentum space, 
this corresponds to the breakdown of the perturbative 
treatment of the /cr-dependence at small- A:t. It is in this 
region that the concept of TMD-factorization, incorpo- 
rating TMDs with intrinsic non-perturbative transverse 
momentum, becomes very important. 

While the ^T-dependence at large bp cannot be calcu- 
lated directly in pQCD, the scale dependence can still 
be handled with the evolution equations (|18II22[) . But a 
prescription is needed for matching the large and small 
6T-behavior. The most common matching procedure was 
developed in Ref. 66] . It replaces b^- in the hard part of 
the calculation by a function. 



b,(b2 



br 



v/1 + b^/b 



2 

max 



(24) 



This definition of b*(bT) is constructed so that it is equal 
to hx when is small, while smoothly approaching an 
upper cutoff 6max when b^ becomes too large. The value 
of 6max is typically chosen to be of order ^ 1 GeV~^ and 
should be thought of as characterizing the boundary of 
the perturbative region of the br — dependence. 

In the calculation of the hard coefficient in Eq. (|23p . 
the appropriate size for the scale ^ is determined by the 
size of bH.(bT). Hence, we define the scale. 



Mb 



Ci 



b^hr) 



(25) 



The parameter Ci is chosen to optimize the perturbation 
expansion. For all our calculations, we will use Ci — 
2e At large br in the final expression for the evolved 
TMD PDF, the effect of the deviation between b^ and 
b* in Ff/p{x, b^r; fJ-', Cf) and K{hT', /i) will be accounted 
for by extra non-perturbative, but universal and scale- 
independent, functions. 

Applying the evolution equations in Eqs. (I18ll22p . using 
the coUinear factorization treatment for small by from 
Eq. (|23[) . and implementing the matching procedure of 
Eq. (I23D allows the TMD PDF to be written with maxi- 
mum perturbative input in terms of evolution from fixed 
starting scales: 
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^ dx ~ 



exp In K{bt;^b)+ / —r 



7i.(5(/i');l)"ln- 



7-7k(5(m')) 



X exp gj/p{x, br) + gKibr) In ■ 



(26) 



This is our master equation for fitting TMD PDFs while incorporating evolution. For a much more detailed explanation 
of the steps summarized above and leading to Eq. ([26)) . we again refer the reader to Ref. [2^, especially chapters 10 
and 13. The steps for evolving are very similar to traditional applications of the CSS formalism, but now they are 
applied to separate, individual TMDs. The scales used in the evolution are chosen to minimize the size of higher 
order corrections in the perturbatively calculable parts. We have labeled three separate factors by "A", "B", and 
"C" to aid in the detailed discussion that will appear in the next section. The C//j(x/a;, 6*; /i^, /ib, K(b„;^b), 
jp{g{fi');l), and ^KigifJ-')) functions are all perturbatively calculable for all by. They are provided to order Ug 
in Appendices. El and IB] On the first line, fj/p{x,iJ,b) is the ordinary integrated PDF from collinear factorization. 
The functions gj/p(x, br) and gKibr) describe the non-perturbative bx-behavior in Ff/p{x, b^; /i; Cf) and K(hT', /i) 
respectively. They are scale- independent and universal. The function gKibr) is notably independent of the species 
of external hadrons. Our definition of the factor gj/pix^br) differs slightly from what is used in [2^ because it has 
absorbed a term equal to gKibr) ln(-y/C^o/a;Mp). This will allow us to choose an arbitrary starting scale (^p^ for the 
evolution in VCF- 

We are ultimately interested in the momentum space TMD which is just the Fourier transform of the coordinate 
space TMD PDF in Eqs. (|12I26I) : 



Ff/pix,'kr;fiXF) 



i2n) 



d^hre'^^^-"^ Ff/pix,hr;t^,CF 



(27) 



Once Eq. (1261) has been parametrized, Ff/pix, k^; ^J■^ Cf) can be determined directly by a numerical Fourier transform. 
Exactly analogous steps hold for the TMD FF. It is related to Kihr', m) by. 



dlnDiz,hr;n,CD) 



Kihr;^l). 



There is also an RG equation analogous to Eq. ()2ip . with anomalous dimension 7D(<7(/i); Cij/m^)' 

dlni)(z,bT;M:C-D) 



dln/i 



7d(5(m);Cd/m )• 



For the small-br region, the collinear factorization treatment of Eq. (fT4|) . analogous to Eq. (|23p . gives 

^ /■! dz 

DH/fiz,hr;fi,(;D) = X! / ^7^5 Cd, /i,5(M))4/j(2; m) + C'IIAqcd^t)") 



(28) 



(29) 



(30) 



The analogue of Eq. ([261) for the TMD FF is 



^ff//(2^,bT;M,Ci 



E 



dz 



Co/fiz/z,b^;^il,lib,gi^ib))dH/]iz,lJ.b) 



X exp { In Kib^ ; ^h) 
^J.b 



dfj/_ 



7-d(5(m'); 1) - In -^^7k(5(m')) 



X exp I gh/jiz, br) + gKibr) In 



V Cd,o 



(31) 



As with the TMD PDF, the perturbative parts of Eq. ((25)) have been calculated to order and are supplied for ref- 
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erence in the appendices. (Note the factor of z^"^"^ that 
appears in the integration measure in Eq. pip as com- 
pared to the factor that appears in Eq. ([26l) : this 
is due to differences in normahzation of the integrated 
PDFs and FFs). The non-perturbative function gKibr) 
is the same in both the TMD PDF and FF. The func- 
tion gh/j{z,bT) describes the non-perturbative large-&T 
behavior that is specific to a fragmentation function for 
parton j and hadron H. The momentum-space TMD FF 
is defined to be 

DH/f{z,zkT;n,CD) = 



d'hr e-^''-^- DH/fiz, br; Co). (32) 



Note that the standard momentum space definition has 
zkr as the transverse momentum argument rather than 

The important resuh of this section is that we now 
have expressions for the evolved TMD quark PDF and 
FF that can be used in Eq. ([2]) , which in turn has a very 
similar structure to the generalized parton model picture 
in Eq. (P). 



VI. IMPLEMENTING EVOLUTION 

We now discuss explicit calculations of evolved momen- 
tum space TMDs. Given some non-perturbative input for 
the large-^T behavior at some fixed scales, we can calcu- 
late the TMDs at different scales by directly calculating 
Eqs. p6H3T|) . We will discuss the TMD PDF and FF 
cases separately. 



A. TMD PDFs 

We first analyze the TMD PDF by discussing each fac- 
tor labeled in Eq. (f26| separately. The first factor (the A- 
factor) matches the TMD PDF to a collinear treatment 
in the small 6t << 1/Aqcd limit. As with standard 
collinear factorization, it involves a hard part, which in 
this case is the coefficient function Cf /j, and a collinear 
factor, which is just the standard integrated PDF. At 
lowest order in a calculation of the coefficient function, 
the A-factor is simply /(a;, /if,). The first factor on the 
second line, the B-factor, is an exponential of quantities 
that can all be calculated perturbatively. They are the 
CS kernel K{b^:;^b) at small br, the anomalous dimen- 
sion 'yp of the TMD PDF, and the anomalous dimension 
7k of the CS kernel. The last factor, the C-factor, im- 
plements the matching between the small and large bp- 
dependence. The function gj/p{x,bT) parametrizes the 
non-perturbative large-&r behavior that is intrinsic to the 
proton, while gKibp) parametrizes the non-perturbative 
large-b^ behavior of K(hp] fi). The function gj/p{x, bp) 
is universal, but in principle depends on the external 



hadron. The function gK{bp) is both universal and in- 
dependent of the species of external hadrons. Note that, 
while the description of the fo^-behavior becomes non- 
perturbative at large bp, there is still perturbatively cal- 
culable evolution for the TMD coming from the B-factor. 

For doing calculations, a choice for the numerical val- 
ues of Ci? and (^p) in Eqs. (|26|3ip is needed. Since 
VC-fC-D ~ Q^i we will treat the PDFs and FFs sym- 
metrically and use -v/Cf = VCd — Q- (In principle, 
slightly different choices may be preferred in specific ap- 
plications, but this will be sufficient for now.) Also, we 
relabel -v/Cf^ 



\/ Cd,o = 2Qo- 
It is instructive to investigate the relationship between 
the parton model expectation and Eq. (|26p . In stan- 
dard collinear factorization for processes integrated over 
transverse momentum, the parton model description of 
the integrated PDF is recovered by dropping all order-as 
contributions to the DGLAP evolution kernel, reproduc- 
ing the Bjorken scaling property of the parton model. In 
collinear factorization, the parton model can be under- 
stood as the zeroth order contribution to the full pQCD 
factorization result. In the TMD PDF case, however, 
if all order-tts or higher contributions to Eq. (j26p are 
dropped, then the TMD PDF becomes: 

Fj/p{x,hp;Cp,fi) -> fj/p{x) 

X exp l^gj/p{x, bp) + gxibp) • (^3) 

Usually, a Gaussian model is used in a partonic de- 
scription of the TMD PDF like Eq. (P). So we write 
gj/p{x,bp) as —gib^/2 and gK{bp) as —g2b^/2. Then 
Eq. ([33)) becomes 



fj/p{x) X exp 



.91 + 52 In 



_Q_ 

2Qo 



bl 



(34) 



This is almost the Gaussian/parton model form of the 
TMD PDF. However, there is still scale dependence com- 
ing from the coefficient of the gxibp) function. This 
difference from the collinear case is due to the fact 
that, while the DGLAP evolution kernels vanish when 
order-cks terms are neglected, the evolution kernel in 
Eg. pJ)) is non- vanishing at zeroth order because of 
the non-perturbative contribution at large bp. TMD- 
factorization therefore differs in a significant qualitative 
way from collinear factorization in that the naive expec- 
tation from the parton model picture is not exactly re- 
covered even in a zeroth order treatment — there is still 
potentially large scale dependence at large bp. This can 
have a large effect on the small- kp scale dependence of 
the TMDs, as already noted in Ref. [13, SJ. In partic- 
ular, if 5i ^ 52, then it can be seen from Eq. ([M]) that 
the TMD PDF becomes extremely sensitive to Q near 
Q ^ 2Qo and at large bp. In the momentum space TMD 
PDF, the evolution corresponds to rapid suppression at 
small kp, of order kp ^ 1 GeV, with increasing Q. The 
effect can be observed in the small-fc^ region of the curves 
in Fig. m 
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Once the A and B and C factors are known, it be- 
comes straightforward to calculate the Fourier transform 
in Eq. (|27|) . Of these, the A- factor is the most cum- 
bersome to deal with because it requires numerical inte- 
grals over X that involve integrated PDFs. The integrated 
PDFs themselves need to be imported from previous fits. 
In our calculations, we obtain the A-factor in Eq. (pS)) 
by using the MSTW PDFs [6Z], along with the MS coef- 
ficient functions calculated in Appendix [X] To facilitate 
future calculations, we have made separate tables for the 
A-factor available for each quark flavor [l[. The B-factor, 
up to order , is straightforward to calculate directly us- 
ing the anomalous dimensions provided in Appendix |B] 

All that is then needed to obtain Eq. (|27|) is a model 
or a fit of the non-perturbative bT-behavior of the C- 
factor. For our calculations, we appeal to currently avail- 
able fits. In principle, fitting the non-perturbative parts, 
(7jyp(a;, 6t) and gKibx), requires knowledge of the com- 
plete (x, b^) plane at different values of Q and for each 
flavor. There have been extensive efforts over the past 
several decades to determine these parameters from ex- 
periments, most commonly from fits to DY processes. 
Currently, the most detailed global fits use the Brock- 
Landry-Nadolsky-Yuan (BLNY) form for the full non- 
perturbative bT-dependence, which leads to a factor in 
the full cross section equal to [39| : 



exp 



91 + 92 In 



_Q_ 

2Qo 



gigs ln(100xAa;B) 



(35) 

The variables xa and xb are the usual momentum frac- 
tion variables of the annihilating quark and anti-quark. 
This almost gives the simple form in Eq. (p4)) . but 
now there is a term in the exponent with explicit x- 
dependence. In the pp cross section, two C-factors ap- 
pear; one with a function 9j/p{x,bT) for the probabil- 
ity of finding a quark in a proton, and the other with a 
function gj^p^x, bx) for finding an anti-quark in an anti- 
proton. Assuming flavor independence, the symmetric 
role of the PDFs in the DY factorization formula allows 
for an immediate identification of the C-factor contribu- 
tion to the TMD PDF in Eq. ([26l) : 



exp 



^In^ 
2 2Qo 



91 



1 



gsMiOx) 



b^ 



(36) 



The fits of Ref. (s^ found gi = 0.21 GeV^, 32 = 
0.68 GeV^ and 53 = -0^6, using Qo = 1.6 GeV us- 
ing data from Refs. [68l - l73| . However, these fits mix 
data for pp, pp and pCu scattering which means that 
it must be assumed that the non-perturbative func- 
tions gj/p{x,bT), g-j/p(x,bT) and gj/cui^jbr) are simi- 
lar. This is not a serious problem at large Q because then 
the Q-behavior comes mainly from the gKibx) function 
which is independent of external hadrons. However, we 
also want our TMD PDF to be valid at smaller Q ^ Qo 
scales, relevant to many SIDIS experiments. Our strat- 
egy then is to match the BLNY fit to the recent scale- 
independent Gaussian fits by Schweitzer, Teckentrup and 



Metz (STM) [335. Using HERMES SIDIS data [zl IzE 
for {x} = 0.09, (g2) w 2.4 GeV^ and z > 0.2 they find 



Ff/p{x,kT) = ff/p{x) X 



exp [—k^/ik^)] 



(37) 



with = (.38 ± 0.06) GeV^. To recover this in our 
fit, we modify the BLNY parametrization in Eq. p6p by 
rewriting it as. 



exp 



^In. 



91 



2Qo 
1 
2 



53 In 10 



XXq 



Xq 



bl 



(38) 



If xo ~ 0.02, then Eq. ((38|) approximately matches the 
STM fit for x = 0.09 and Q = y/2A GeV, but reduces 
to the BLNY fit at larger Q and smaller x. We note 
that the x and b-p dependence does not quite factorize 
in these TMD fits. Indeed, the form of the gj^p{x, bx) is 
not required by the formalism to factorize into separate 
X and by dependence. 

We now have a fit that includes the scale dependence of 
the QCD evolution in Eq. (^5)) . and whose ^T-dependence 
matches two previously performed fits for different re- 
gions of kinematics. For illustration, we have plotted 
in Fig. [U the TMD PDF of the up-quark for the small, 
medium and large values oi Q = -\/2.4, 5, and 91.19 GeV 
and with x — 0.09. We have made the plot run over 
a range from /cy = to 6 GeV, typical for studies of 
TMD-functions. (Recall, however, that without the Y- 
term of Eq. © the TMD PDF by itself only has a simple 
interpretation for kx << Q-) Comparing the curves, it 
is clear that the evolution in Q is a large effect, lead- 
ing to more than an order of magnitude of suppression 
at small kx, and a broad tail at larger kx- Numerical 
computations that produce plots like Fig. [1] are avail- 
able at Ref. [ij. The large bx cutoff, 6„iax, should 
be small enough to exclude the non-perturbative large 
^T-regime from the perturbative part of the TMD PDF, 
but it is otherwise arbitrary. In fits, different choices 
of braax Can lead to very different values for the non- 
perturbative parameters because changing the size of 
^max effectively reshuffles contributions to the TMD PDF 



between the different factors in Eq. ([26|) . It turns out 
that the b^ax = 0.5 GeV'^ value from the BLNY fit is 
rather small. In other words, this choice of 6niax restricts 
the perturbative part of the calculation to a range in bx 
that is significantly smaller than the range where per- 
turbative methods are still reasonable. The analysis in 
Ref. [76] has found that 6max = 1-5 GeV~^ is preferred, 
and the parameters for the BLNY fit in Eq. (|35p become 
gi = 0.201 GeV2, 52 = 0.184 GeV^ and 33 = -0.129. 
By using the newer parameters from Ref. [76|, we can 
again construct a TMD PDF parametrization from the 
BLNY form that matches the STM fit at small Q by us- 
ing Eq. ([551) . With the newer parameters we find that 
Xq — 0.009 is needed to fit to the STM parametrization 
at small kx- 
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Up Quark TMD PDF, x = .09 
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FIG. 1: The up quark TMD PDF for Q = ^/2l^, 5.0 and 91.19 GeV and x = 0.09. The upper plot shows the result of using the 
BLNY fit in Eq. (|38j with b^nax = 0.5 GeV'^ while the lower panel shows the BLNY fit obtained with bmax = 1-5 GeV"\ The 
solid maroon, dashed blue, and red dot-dashed curves are for Q — \/2.4, 5.0 and 91.19 GeV respectively (see online version for 
color). 



Up Quark TMD PDF, x = .09, Q = 91.19 GeV 




FIG. 2: Comparing the shape of the TMD PDF within various approximations. The solid red curves are the same as the 
Q — 91.19 GeV curves in Fig. [T] The dashed blue curve is the result of setting the A-factor in Eq. (|26[) equal to f{x,iJ,b), and 
the dash-dotted maroon curve is obtained by setting the B-factor in Eq. (|26p equal to 1. (See online version for color.) 
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One reason that we prefer the smaller 6max for the 
present paper is that our present analysis includes only 
the order-a^ contributions to the perturbatively calcu- 
lable parts, so it is important that higher order con- 
tributions are small. In practice, higher order contri- 
butions can have a large effect, especially at small and 
intermediate kr- We estimate the size of the theoreti- 
cal error in our analysis by redoing the calculation for 
the parametrization with the larger value of &jnax = 
1.5 GeV~^, and using the parameters of Ref. [7a]- The 
result is the lower panel in Fig. [TJ By comparing the up- 
per and lower plots, it can be seen that the curves differ 
by a maximum of about a factor of two for the large Q 
curve. (By running the calculation for different values of 
X, we have verified that this is generally true for x be- 
tween about 0.01 and 0.2.) The largest difference is for 
the Q = 5 GeV curves. This is because for Q = 5 GeV 
evolution effects become significant, but the different val- 
ues of xq needed for the curves to match at Qo lead to a 
significant difference in fc-p dependence. 

In the future, it will be possible to decrease the theo- 
retical uncertainty in the TMD PDF parametrization by 
including higher orders in the perturbative parts of the 
calculation, and by using improved fits based on newly 
available data. 

It is also instructive to investigate the effect of the sep- 
arate A, B and C factors in Eq. (|26l) . In Fig. [5J we have 
again plotted the TMD PDF for Q = 91.19 GeV. In ad- 
dition, we show the effect of replacing the A factor by 
simply the lowest order, unevolving result, /(s:, /ib^^x) 
(blue dashed curve), and the effect of replacing the B 
factor by one (maroon dash-dotted curve). (The color 
version is online.) The dashed curves show that simply 
using /(a;, A*6max) instead of the full A factor is typically 
a very good approximation in the small-fc^ limit of the 
TMD PDF. This significantly simplifies the calculation 
of the TMD PDF in cases where the very small-fc^ re- 
gion is the main contribution of interest. (However, it 
should be re-emphasized that the full A factor is needed 
for a complete description of the cross section over all 
qr up to order Q.) Neglecting the B factor introduces 
a substantial error at small kx and completely removes 
the large-fcr tail. We have compared the calculations for 
the two different 6max values, 0.5 GeV^^ and 1.5 GeV~^, 
in the upper and lower plots. As could be expected, 
the effect of setting _B = 1 is substantially larger in the 
&max = 1.5 GeV~^ case where the role of higher orders is 
more important. 

It is common in TMD studies to use Gaussian 
parametrizations like Eq. (I37p . An attractive feature of 
such an approach is that the TMD has a simple and well- 
defined integral over all kx, and the standard integrated 
PDF is obtained simply by integrating the TMD PDF. 
Moreover, the Gaussian form makes the calculation of 
weighted structure functions simple. Therefore, it is use- 
ful to investigate how well a Gaussian form describes the 
shape of the TMD, and over what range of transverse mo- 
mentum. (However, we recall that the actual relationship 



between integrated and TMD PDFs is more complicated, 
as already evidenced by the broad tail in Fig.[2j) To give 
an example of such a comparison, and to study the ef- 
fect of the tail in fits of the TMD PDF at large Q, we 
have replottcd the Q = 91.19 GeV curves from Fig. [T]in 
Fig. [3l but now we include Gaussian fits. From the plots 
it can be seen that at large-Q the Gaussian shape contin- 
ues to do a reasonable job of describing the very small fcr 
behavior (less than a few GeV), but since it completely 
neglects the tail at large fcr tail it underestimates the size 
of the typical kx- 

To investigate the role of the tail, we calculated the 
integral of the TMD PDF over kx, weighted by k, 



kl= I d^kTk^F{x,kT), 



(39) 



and we compared the result of using the Gaussian fit with 
the result obtained by numerically integrating the orig- 
inal Q = 91.19 GeV curve (the solid curve in Fig. [3|). 
For the original curve, Eq. (|39|) is quite ill-defined be- 
cause of the large-A:T tail. In addition, the contribution 
from very large kx is outside the region where a TMD- 
factorization description alone is valid, and the y-term 
becomes important. Nevertheless, it is possible to get 
a sense of the effect of the tail on typical values of kx 
by integrating up to a cutoff that is large, but still sig- 
nificantly less than Q. For the TMD PDFs in Fig. |3 
we choose an upper cutoff of 20 GeV. For the Gaussian 



fit, we find 



'fc2 = 6 GeV for 



= .5 GeV-i and 



fc2 = 4 GeV for 6„ 



^max = 1.5 GeV ^. (The difference is 
due to the slightly different ranges in kx where a Gaus- 
sian is a good fit.) For the original curve, with the upper 



cutoff on kx of 20 GeV, we find 

= 0.5 GeV-i and b^,., = 1.5 GeV'i 



k^ = 15 GeV for both 



Hence, for 



both values of b„ 



the tail leads to at least a factor of 



two increase in the typical kx- 



B. TMD FFs 

The non-perturbative input for the FFs is much less 
constrained by existing analyses. However, the function 
gK{bx) in Eq. ^ for the TMD PDF is the same func- 
tion that appears in Eq. ^ for the TMD FF. There- 
fore, given a fit for the TMD FF at a particular scale, 
one can use the same gxibx), along with the anoma- 
lous dimensions and coefficient functions calculated in 
Appendices lAIBl to estimate the evolution to different 
scales. For the starting scale, we again appeal to the fit 
of Ref. which uses a Gaussian form, 

DH/fiz, Kt) = dH/f{z) X J^ , ^ (40) 

where Kx is the hadron transverse momentum in the 
photon rest frame. Again fitting the HERMES data. 
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Up Quark TMD PDF, x = .09, Q = 91.19 GeV 





FIG. 4; TMD FF for charged pions from a hadronizing up-quark. The upper plot is for g2 = .68 GeV^ and the lower plot is 
for g2 = .184 GeV^ In each case, TMD FF matches to the STM fit at Q = V^A GeV. (See online version for color.) 
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for SIDIS in the kinematical range, (x) = 0.09, (Q^) w 
2.4 GeV^ and z > 0.2, they find that (K^) = (.16 ± 
0.01) GeV^. One can write the transverse components of 
the photon-frame hadron momentum Kt in terms of the 
transverse components of the hadron-frame parton mo- 
mentum kr as Kt = —zkr- The analogue of Eq. (p4|) 
for the FF has an extra so that when aU order as 
corrections are dropped the FF reads 



DH/f{z,bT-XF,lJ') 

X exp 



1 



' 1 2 1 Q 
91+92Z In:^ 



Equating this to the inverse Fourier transform of 
Eq. PU)) . we identify the factor in brackets as 



9i + 92Z 



2 In 



V2.4GeV 
2Qo 



0.08 GeV. 



From this relation we can extract a value for g[. The 
factor multiplying —b^ in Eq. (|4T1) can then be identified 
with the non-perturbative exponential factor in Eq. (1311) . 
Using Refs. [77h80 | for the integrated FFs, we can then 
calculate the TMD FF using Eq. We have repeated 
the analysis of the TMD PDF for a TMD FF of a charged 
pion fragmenting from an up quark. Fig. 2] shows the 
TMD FF for different energy scales, Q = ^/2A , 5 and 
91.19 GeV. By comparing different energy scales, one can 
immediately see the effect of including perturbative evo- 
lution in the definitions of the TMD FFs from the high 
kr tails the TMD FFs acquire. We have also repeated 
the analysis of evaluating the TMD FF for different val- 
ues of 5max = 0.5 and 1.5 GeV~^ and we find a similar 
error estimate as in the case of the TMD PDF. The com- 
parison is shown again in Fig. S) Note that in Fig. U] we 
have plotted the TMD FF as a function of the hadron 
transverse momentum Kt rather than parton transverse 
momentum kx- 

We also investigated how well a Gaussian function fits 
the perturbatively evolved TMD FF. As with the TMD 
PDF, the Gaussian fit does not adequately capture the ef- 
fects of perturbative evolution for the TMD FF. The con- 
tribution of the fcr-tail is smaller in the case of the TMD 
FF. This can be understood by comparing the fc^ depen- 
dence of a TMD PDF with a TMD FF. The TMD FF is 
less broad in fc^ than a TMD PDF and therefore drops 
faster with a smaller fc^ tail. To quantify this we have 
once more calculated a typical fc^ using Eq. (15^ both for 
the Gaussian fit and the actual TMD FF. For 6niax = 0.5 



GeV ^ we find that for the Gaussian fit 



1.74 



GeV while for the actual TMD FF J k^ = 2.15 GeV 
which gives a relative difference of 23.5%. For the case 
of fcmax = 1-5 GeV~^ the values are -y/^ = 1.06 GeV 

for the Gaussian fit and = 1.85 GeV for the actual 
TMD FF with a larger relative difference of 73.5%. 



VII. DISCUSSION AND CONCLUSIONS 

Factorization theorems provide the bridge between 
abstract field theoretical concepts and phenomenology, 
and are responsible for giving pQCD its great predic- 
tive power. The parton distribution and fragmentation 
functions, which arise naturally from the factorization 
derivations, play a central role in relating formal pQCD 
to parton model concepts. A precise understanding of 
the definitions, evolution, and universality properties of 
these correlation functions is what enables calculations 
in pQCD to make accurate first principles predictions. 

While the standard formalism of collinear factoriza- 
tion has proven extremely useful for sufficiently inclu- 
sive processes, the more sophisticated formalism of TMD- 
factorization is needed for processes in which the intrinsic 
transverse momentum of the partons becomes important. 
As has already been widely discussed, there are a number 
of technical and conceptual subtleties involved in arriv- 
ing at good definitions for the TMDs that are consistent 
with factorization. These issues include the need to reg- 
ulate and deal with rapidity divergences and achieve a 
cancellation of spurious Wilson line self-energies. The 
subtleties involved in defining TMDs have largely been 
clarified and resolved in Ref . :26] , which provides defini- 
tions that are consistent with the requirements of factor- 
ization, and demonstrates the relationship with the usual 
CSS formalism. 

While considerable effort has been devoted to imple- 
menting CSS evolution in unpolarized scattering, the re- 
sulting parametrizations are often not framed in the lan- 
guage of TMD PDFs. By contrast, for polarization de- 
pendent TMDs, there has been very little work done in 
implementing evolution in the parametrization of exper- 
imental data. Up to this point, these functions have only 
been probed over a very narrow range of scales so that 
evolution has not been a major issue. However, for future 
progress in understanding the role of quark and gluon 
degrees of freedom in hadronic structure, it will be im- 
portant to remedy this situation. Ideally, there should 
be collections of tabulated fits to the TMDs that incor- 
porate evolution, and which can be directly related to 
the field-theoretic definitions of the correlation functions, 
analogous to what has already existed for some time in 
collinear factorization. 

We have started this process by recasting previously 
performed fits [H,!!!!!!] of unpolarized TMD PDFs and 
TMD FFs in terms of the TMD definitions of Ref. [H]. 
This provides a much clearer connection between the for- 
malism of evolution and generalized parton model ap- 
proaches, and provides practical TMD parametrizations 
that can be used directly in TMD calculations. We have 
also completed the derivation of the lowest order anoma- 
lous dimensions and coefficient functions for the TMD 
PDF. At our website [l[, we have supplied tables and 
interpolation routines for the parts of the quark TMD 
PDFs and FFs that can be described using collinear fac- 
torization (the "A-factors" in Eqs. \ST\i ] for each fla- 
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vor of quark, as well as sample calculations that give plots 
like Fig. □ 

We have confirmed the important observation that evo- 
lution has a strong quantitative effect on the TMDs and 
therefore should be included in future phenomenological 
applications of TMD factorization, particularly given the 
range of energy scales that are set to be probed in the 
future. Another reason to have reliable fits of TMDs, 
including evolution, is that it opens the possibility to 
identify instances of factorization breaking effects of the 
type discussed in Ref. 45-48]. Recognizing factorization 
breaking effects will be an important next step in expand- 
ing our understanding of pQCD phenomenology. Even 
in unpolarized scattering, there is a possibility to use 
parametrizations like those presented in this paper to test 
the factorization hypothesis. Recent RHIC data f57l.[58|. 
for example, may be useful for such an analysis. 

Nevertheless, much work remains to be done. The the- 
oretical uncertainty in the TMD fits can be reduced by 
including higher orders in the calculations of the anoma- 
lous dimensions, the K kernel, and the collinear coeffi- 
cient functions. Moreover, as new data from experiments 
like those taking place at the LHC, RHIC, JLab, and a 
possible electron-ion collider are made available and ana- 
lyzed, it will be possible to obtain improved fits. Already, 
there are data from ATLAS [8l| which can potentially 
help to improve the quality of fits for the unpolarized 
TMDs. TMD effects can also be studied in an e+e" 
collider as recently discussed in Ref. 41|. A number 



of theoretical issues with the evolution formalism itself 
also remain unsettled. For instance, the precise form of 
the matching function for between perturbative and non- 
perturbative transverse momentum regimes in Eq. (|24l) is 
somewhat arbitrary and better prescriptions may be pos- 
sible. Along similar lines, a truly optimal choice of hmax 
may be different from the values we have used here. One 
possiblility may be to formulate the evolution directly in 
momentum space. 

One of the most important next steps is to extend the 
analysis presented in this paper to the Sivers and Boer- 
Mulders functions, which are needed for clarifying the 
spin structure of hadrons. Efforts to address polariza- 
tion dependent situations can utilize existing fixed-scale 
fits (such as [H-IH, [HI). In such cases, a careful treat- 
ment of the matching between large and small transverse 
momenta will also be important 's^]. Furthermore, it 
will be important to establish the relationship be evolved 
TMDs and the evolution of weighted functions such as 
those treated in Ref. [43| . 

Fits of the gluon TMD PDF that include CSS evo- 
lution will be also be needed, especially in tests of fac- 
torization. See [S^ [83| for recent work related to evo- 
lution and gluon resummation in the context of gluon 
PDFs. In addition, recent calculations in Ref. [85| have 
shown how to probe linearly polarized gluons in heavy 
quark production, and the universality properties of the 
gluon TMD PDF have been clarified in Ref. For 
processes that probe the gluon TMDs, some important 



details of the TMD-factorization theorems have yet to 
be completely understood. The issue of so-called "super- 
leading regions" in the factorization theorems that use 
the gluon distribution [s?] still needs to be clarified in 
the TMD case. Furthermore, for processes that involve 
several final state hadrons, such as e+p Hi + H2 + X , 
the separation of the soft factor into universal square 
root factors as in the definitions in Eqs. and 
is not straightforward. Following a naive analysis like in 
Sect. IIVCI seems to suggest that extra soft factors are 
needed, and that a more complicated factorization struc- 
ture is required. 

These are all issues we intend to pursue in a continua- 
tion of the TMD project. 
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Appendix A: Coefficient EYinctions 



present the steps for cal- 
coefficient functions (the 
functions in the A-factor 



In this appendix, we 
culating the collinear 
Cf/j{x/x,K;nl,fii„g{Hb)) 
of Eq. ^W)- We first briefly review the steps, presented 
in Ref. 1^, for the case of the TMD FFs. Then we 
explain the extension to the analogous case for the TMD 
PDFs. 

In perturbation theory, the FF in Eq. (IT^ itself obeys 
a collinear factorization theorem [ij, valid for small b^. 
The collinear part is just the standard integrated FF. 
Writing the factorization a.s D = d (E) C , we have to first 
order 



dW ^ SO] ^ (5(1] _^ ^[1] ^ (^[0] 



(Al) 



The superscripts label the order in perturbation theory. 
Using d!jj-,{z) — Sjj'S{z — 1) for the lowest order inte- 
grated FF, one finds 



,2-2£ 



(A2) 



for the first order FF coefficient function. To get the 
collinear coefficient function, all that is needed then is to 
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calculate the first order expression for the unintegrated 
FF D^j^f (2 , by) and for the integrated case dj/f{z). Since 

C^^Jj{z,hT) is independent of the species of external 
hadron, the calculation can be done for the special case 
of a quark hadronizing to a gluon. 



The order 0{g'^) diagram is shown in Fig. [5] There 
is no leading contribution from the soft region and hence 
no need to subtract soft factor contributions. Calculating 
the first order TMD FF gives 



(2^) 



4-2e . 



87r3 



1 + (1 - z)2 - 



(A3) 



The corresponding integrated FF is calculated in nearly the same way, except that hr is set to zero and the 1/z factor 
in the definition (|13p is changed to z^^~^ in the integrated case. Also, an MS counterterm is needed to remove the 
resulting UV divergence. The result is 



d'^j (z) = 



87r3 



l + (l-z) 



2 



g'gF(47r)^ 
87r2r(l - e)e 



l + (l-z)2 



(A4) 



Performing the integrals and putting all the terms together in Eq. (|A2[) gives the coUinear coefficient function 

agCp 



Cg/f{z,hT]iJ.; Colli ) 



2[l + (l-zf] 



(A5) 



For the TMD FF of a quark hadronizing to a quark, the diagrams are shown in Fig. [HI along with the soft subtraction 
factors in Fig. [7] that are needed according to the definition in Eq. p^ . Apart from the need to include soft factor 
contributions, the steps are analogous to those that led to Eq. (|A5p The integrated FF is found again using Fig. [6] 
along with MS counterterms. The result is 



1 



27r 



In 



2z 



7e 



1 - z 



1 



1 



+S{l-z) 



2 [In {b■'^^l^) - 2(ln2 - 7e)] - [Hbrl^^) " 2(ln2 - 7e)] In ^ 



+ Oiai). (A6) 



Again, the details for the above FF calculation can be found in Ref. [26|. 

The steps for calculating the order-C(52) contributions to the collinear coefficient functions for the TMD PDFs 
are analogous to the steps used in the TMD FF case, with some minor changes. We provide them here, presented 
for the first time in the context of TMDs. The results should be equivalent to calculations already done in the CSS 
formalism, up to possible changes in scheme. For the TMD PDF, the analogue of Eq. (jA2p is 



(7[)i,(.,b,) = i?j;i(x,b,)-/|;i(x). 



(A7) 



The difference in factors of longitudinal momentum fraction in the second term comes from the different normalization 
in Eq. ([50]) as compared with (1^51) . Again, we may perform the calculation for on-shell external partons. Using the 
diagram in Fig. [SI for the TMD PDF of a quark inside a gluon one has 



2,,2e 



(x hr) = 



(2^)4~2e(l_£) 



dk-cP~^"kT €-'^-^-^^2715 {{p ~ kf) 



l SpoiTr(7+|^^,(^~;^)0 
4 



(fc2)2 



8^2r(i - e) J, 



1 - 



2a;(l - x) 



(AS) 



where is the polarization vector for the initial state gluon and we sum and average over all possible polarizations. The 
integrated PDF is found again by setting = in the above equation and adding the appropriate MS counterterm 
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FIG. 5: One-loop diagram used contributing to the TMD FF and the integrated FF of a quark fragmenting to a gluon. 

/ / 

k 




FIG. 6: One-loop diagrams contributing to the TMD FF and the integrated FF of a quark fragmenting into a quark. 



for the resulting UV divergence. This gives 



[1] _ .9^(4V)-T, 



2x{l - x) 
1-e 



87r2r(l - e)e 



l-2x{l-x) 



(A9) 



Using Eas. (jA7P - (IA9p and evaluating the kx integrals, gives the TMD PDF coUinear coefficient function for finding a 
quark of flavor j' in a gluon at order as, 



Cj7g(a;,bT;Ai;CF/Ai^) 



2tt 



2 [1 - 2x{l - x)] 



7E 



+ 2x{l-x) +0{ai) 



(AlO) 



Finally using diagrams in Fig. [9] together with the soft subtraction terms in Fig. [7] for the TMD PDF for finding a 
quark of flavor j' in a quark of fiavor j and again the diagrams in Fig. [Altogether with the MS UV counterterms for 



Vn 



Vn 




FIG. 7: One-loop diagrams for the soft-factor contributions of Eq. (I31|) . Hermitian conjugate graphs are also needed but are 
not shown. 
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k-p 




FIG. 8: One-loop diagram contributing to the TMD PDF and the integrated PDF for a quark inside a gluon. 




FIG. 9: One-loop diagrams contributing to the TMD PDF and the integrated PDF of a quark inside a quark. Hermitian 
conjugate graphs are also needed but are note shown. 



the integrated quark PDF we find to order as 



27r 
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l-x 
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'C 
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Appendix B: Anomalous Dimensions 

All calculations of anomalous dimensions defined in 
Eqs. (|20 |) . pip and ([29l) use dimensional regularization 
with the MS scheme. The anomalous dimension of the 
quark TMD PDF up to order Us is, 

7F(M;CF/M') = «.^0-ln(^)) +0(a^). (Bl) 

At order as , the quark TMD FF anomalous dimension is 
the same as for the TMD PDF. We note that these results 
are consistent with what is found in, e.g., Ref. (6l| using 



different methods. 

The CS kernel, in Eq. ([T51) . up to order as hi b^-space 

is, 

k{pL, hr) - [M/i'6|) - ln4 + 27e] + 0{al). 

(B2) 

The anomalous dimension of K (see Eq. (|20p ) is up to 
order as. 



7K(M)=2^+0(af). 

TT 
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